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Abstract. Neutron star configurations are considered as thermodynamical systems 
for which a phase diagram in the angular velocity (O) - baryon number (N) plane is 
obtained with a dividing line iV cr it(J2) for quark core configurations. Trajectories of 
neutron star evolution in this diagram are studied for different scenarios defined by 
the external torque acting on the star due to radiation and/or mass accretion. They 
show a characteristic change in the rotational kinematics when the star enters the 
quark core regime. For isolated pulsars the braking index signal for deconfinement 
has been studied in its dependence on the mass of the star. Model calculations of the 
spin evolution of accreting low-mass X-ray binaries in the phase diagram have been 
performed for different values of the initial magnetic field, its decay time as well as 
initial mass and mass accretion rate. Population clustering of these objects at the line 
iV cr it(i?) in the phase diagram is suggested as an observable signal for the deconfinement 
phase transition if it exists for spinnning and accreting neutron stars. 

1 Neutron stars as thermodynamical systems 

Quantum Chromodynamics (QCD) as the fundamental theory for strongly inter- 
acting matter predicts a deconfined state of quarks and gluons under conditions 
of sufficiently high temperatures and/or densities which occur, e.g., in heavy-ion 
collisions, a few microseconds after the Big Bang or in the cores of pulsars. The 
unambiguous detection of the phase transition from hadronic to quark matter 
(or vice-versa) has been a challenge to particle and astrophysics over the past 
two decades While the diagnostics of a phase transition in experiments 

with heavy-ion beams faces the problems of strong noncquilibrium and finite 
size, the dense matter in a compact star forms a macroscopic system in thermal 
and chemical equilibrium for which signals of a phase transition shall be more 
pronounced. 

Such signals have been suggested in the form of characteristic changes of 
observables such as the surface temperature ||, brightness Q, pulse timing || 
and rotational mode instabilities || during the evolution of the compact object. 
In particular the pulse timing signal has attracted much interest since it is due 
to changes in the kinematics of rotation. Thus it could be used not only to 
detect the occurrence but also to determine the size of the quark core from 
the magnitude of the braking index deviation from the magnetic dipole value 
. Besides of the isolated pulsars, one can consider also the accreting compact 
stars in low-mass X-ray binaries (LMXBs) as objects from which we can expect 
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signals of a deconfinement transition in their interior ■ The observation of 

quasiperiodic brightness oscillations (QPOs) for some LMXBs has lead to 
very stringent constraints for masses and radii which according to p"^ , |l3| 
could even favour strange quark matter interiors over hadronic ones for these 
objects. Due to the mass accretion flow these systems are candidates for the 
formation of the most massive compact stars from which we expect to observe 
signals of the transition to either quark core stars, to a third family of stars 
[|14| or to black holes. Each compact star configuratrion can be identified with a 
thermodynamical system characterized by the total baryon number (or mass), 
temperature, spin frequency and magnetic field as thermodynamical variables. 

Since the evolutionary processes for the compact objects accompanying the 
structural changes are slow enough we will consider here the case of rigid rotation 
only and restrict ourselves to the degenerate systems at T = 0. The magnetic 
and thermal evolution of the neutron stars we will consider as decoupled from 
the mechanical evolution. 

In this approximation we can introduce a classification of isolated and accret- 
ing compact stars in the plane of their angular frequency fl and mass (baryon 
number N) which we will call phase diagram, see Fig. [l]. 

Each point in this phase diagram corresponds to a mechanical state of a 
neutron star. Mechanical equilibrium of thermal pressure with gravitational and 
centrifugal forces leads to a stationary distribution of matter inside the config- 
uration. This distribution is determined by the central density and rotational 
frequency. Thermal and chemical equilibrium are described by an equation of 
state which determines the structure and composition of the compact star con- 
figuration. 

The requirements of stability restrict the region of stability of the quasis- 
tationary rotating star configurations. From the right hand side of the phase 
diagram there is a line which separates a region of stars which are collapsed to 
black holes (BH) from the stable ones. At high baryon numbers beyond this line 
gravitational forces dominate over pressure and centrifugal forces of matter. 

From the top of the diagram the region of stable star configurations is sepa- 
rated by the Keplerian frequencies Qk (N) from that where mass shedding under 
the centrifugal forces does not allow stationary rotating objects. 

Inside the region of stationary rotators one can distinguish two types of 
stationary rotating stars: those with quark matter cores (QCSs) and hadronic 
stars (HS) The region of the QCSs is expected to be located at the bottom 
right of the phase diagram in the region of the most massive and slowly rotating 
compact stars. The critical line N CT it(f2) which separates the QCS region from 
that of HSs is correlated with the local maxima of the moment of inertia with 
respect to changes of the baryon number at given angular velocity fl due to 
the change of the internal structure of the compact object at the deconfinement 
transition . 

Using the analogy with the phase diagram of a conventional thermodynamical 
system we can consider trajectories in this diagram as quasistationary evolution- 
ary processes. 
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Fig. 1. Sketch of a phase diagram for spinning neutron stars with deconfinement phase 
transition. 

Since these processes for neutron stars have characteristic time scales (« 
10 6 -^10 8 yr) much longer than those of typical observers it is almost impossible to 
trace these trajectories directly. A possible strategy for measuring evolutionary 
tracks will be to perform a statistics of the population of the different regions in 
the phase diagram. 

Our aim is to investigate the conditions under which the passage of the phase 
border leads to observable significant clustering in the populations. 

We will provide criteria under which a particular astrophysical scenario with 
spin evolution could be qualified to signal a deconfinement transition. 

2 Phase diagram for a stationary rotating star model 
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2.1 Equation of state with deconfinement transition 

Since our focus is on the elucidation of qualitative features of signals from the 
deconfinement transition in the pulsar timing we will use a generic form of an 
equation of state (EoS) with such a transition Q which is not excluded by the 
mass and radius constraints derived from QPOs. In our case as well as in most 
of the approaches to quark deconfinement in neutron star matter a standard 
two-phase description of the equation of state is applied where the hadronic 
phase and the quark matter phase are modelled separately. The ambiguity in 
the choice of the bag constant for the description of the quark matter phase can 
be removed by a derivation of this quantity |p^| , ^7| from a dynamical confining 



approach 18 . The resulting EoS is obtained by imposing Gibbs' conditions for 
phase equilibrium with the constraints of globally conserved baryon number and 
electric charge f^f^H). 

Since nuclear forces are not fundamental, our knowledge about the equation 
of state for nuclear and neutron star matter at high densities is not very precise. 
There is not yet a description of the equation of state for strongly interacting 
matter on the fundamental level in terms of quark and gluon degrees of freedom 
where nucleons and mesons appear as composite structures. It is one of the 
goals of nuclear astrophysics and of neutron star physics to use pulsars and 
other compact objects as laboratories for studies of the nuclear forces and the 
phases of nuclear matter . The existence of several models for the equation 

of state for dense stellar matter allows a variability in the phase structure of 
neutron stars |^2|,^3), therefore the phase diagram (Fig. |l|) remains robust only 
qualitativly. 

For the detailed introduction of the phase diagram and a quantitative analysis 
using thermodynamical methods, we will employ a particular EoS model [Q 
which is characterized by a relatively hard hadronic matter part. 



2.2 Configurations of rotating stars 

In our model calculations we assume quasistationary evolution with negligible 
convection and without differential rotation which is justified when both the mass 
load onto the star and the transfer of the angular momentum are sufficiently slow 
processes. 

For our treatment of rotation within general relativity we employ a per- 
turbation expansion with respect to the ratio of the rotational and gravitational 
energies for the homogeneous Newtonian spherical rotator with the mass density 
p(0) equal to the central density, E rot /E gYav = (Sl/H) 2 , where Sr = 47rG j o(0). 
This ratio is a small parameter, less than one up to the mass shedding limit 0. 

The general form of the expansion allows us to describe the metric coefficients 
and the distributions of pressure, energy density and hydrodynamical enthalpy 
in the following form 



X{r, 0; Q) = X^\r) + (f2/f2) 2 X^{r, 6) + 0{Q A ) , 



(1) 
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where X stands for one of the above mentioned quantities ffl] ■ The series expan- 
sion allows one to transform the Einstein equations into a coupled set of differ- 
ential equations for the coefficient functions defined in ([!]), which can be solved 
by recursion. The static solutions, obeying the Oppenheimer-Volkoff equations, 
are contained in this expansion for the case J? = when only the functions with 
superscript (0) remain. The other terms are corrections due to the rotation. We 
truncate higher order terms ~ 0(Q A ) in this expansion and neglect the change 
of the frame dragging frequency, which appears at 0(J? 3 ). For a more detailed 
description of the method and analytic results in the integral representation of 
the moment of inertia we refer to JtJ and to works of Hartle and Thorne [p4|^5|] , 
Sedrakian and Chubarian |26| , p7| . 
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Fig. 2. Phase structure of rotating hybrid stars in equatorial direction in dependence 
of the angular velocity O for stars with different total baryon number: Nb/Nq — 
1.3,1.55,1.8. 

In Fig. H we show the critical regions of the phase transition in the inner 
structure of the star configuration as well as the equatorial and polar radii in 
the plane of angular velocity fl versus distance from the center of the star . It 
is obvious that with the increase of the angular velocity the star is deformating 
its shape. The maximal excentricities of the configurations with Nb =1.3 Nq, 
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N B = 1.55 N Q and N B = 1.8 iV Q are e(O max ) = 0.7603, e(J? max ) = 0.7655 and 
e(fi max ) = 0.7659, respectively. Due to the changes of the central density the 
quark core could disappear above a critical angular velocity. 

It is the aim of the present paper to investigate the conditions for a verifi- 
cation of the existence of the critical line N CT i t (f2) by observation. We will show 
evidence that in principle such a measurement is possible since this deconfine- 
ment transition line corresponds to a maximum of the moment of inertia, which 
is the key quantity for the rotational behavior of compact objects, see Fig. 0. 

In the case of rigid rotation the moment of inertia is defined by 

i(n,N) = j{q,n)/q , (2) 

where the angular momentum J{ft, N) of the star can be expressed in invariant 
form as 

J(f2,N) = J Tl^-gdV , (3) 

with being the nondiagonal element of the energy momentum tensor, s/—gdV 
the invariant volume and g = det ||fl^i>|| the determinant of the metric tensor. 
We assume that the superdense compact object rotates stationary as a rigid 
body, so that for a given time-interval both the angular velocity as well as the 
baryon number can be considered as global parameters of the theory. The result 
of our calculations for the moment of inertia (||) can be cast into the form 

/ = I(0)+^4I Q , (4) 

a 

where 1^ is the moment of inertia of the static configuration with the same 
central density and AI a stands for contributions to the moment of inertia from 
different rotational effects which arc labeled by a: matter redistribution, shape 
deformation, and changes in the centrifugal forces and the gravitational field Q . 

In Fig. U we show the resulting phase diagram for compact star configura- 
tions which exhibits four regions: (i) the region above the maximum frequency 
fift(N) where no stationary rotating configurations are found, (ii) the region of 
black holes for baryon numbers exceeding the maximum value iV max (,f2), and the 
region of stable compact stars which is subdivided by the critical line iV cr it(i7) 
into a region of (iii) quark core stars and another one of (iv) hadronic stars, re- 
spectively. The numerical values for the critical lines are model dependent. For 
this particular model EoS due to the hardness of the hadronic branch (linear 
Walecka model j2(J ) there is a maximum value of the baryonic mass on the crit- 
ical line N CT it(Qk) — l-8iV©, such that for stars more massive than that one all 
stable rotating configurations have to have a quark core. This property can be 
seen from the dependence of the phase structure of the star on angular velocity 
in Fig. |[ For the whole interval of possible frequencies in the case of N = 1.8Nq 
the quark core radius remains approximately unchanged: R cor e ~ 7 km. 
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Fig. 3. Phase diagram for configurations of rotating compact objects in the plane of 
angular velocity Q and mass (baryon number N). Contour lines show the values of 
the moment of inertia in Mokm 2 . The line N CI \t(Q) which separates hadronic from 
quark core stars corresponds the set of configurations with a central density equal to 
the critical density for the occurence of a pure quark matter phase. 



3 Evolution scenarios with phase transitions 

We want to explain why the occurence of a sharply peaked maximum for the 
moment of inertia in the fl — N plane entails observational consequences for the 
angular velocity evolution of rotating compact objects. The basic formula which 
governs the rotational dynamics is 

a K(N,n) 

i{N,n) + n{di{N,n)/dQ) N ' {) 

where K = Ki nt +K oxt is the net torque acting on the star due to internal and ex- 
ternal forces. The internal torque is given by K[ nt (N, fl) = —ON (dI(N, f2)/dN) a 
the external one can be subdivided into an accretion and a radiation term 
K eyA = K acc + K Tad . The first one is due to all processes which change the 
baryon number, K acc — N dJ/dN and the second one contains all processes 
which do not. For the example of magnetic dipole and/or gravitational wave 
radiation it can be described by a power law i^ ra d = f3& n , see p8|,p9[. 
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3.1 Spin-down scenario for isolated pulsars 

The simple case of the spindown evolution of isolated (non-accreting, N = 0) 
pulsars due to magnetic dipole radiation would be described by vertical lines 
in Fig. |l|. These objects can undergo a deconfinement transition if the baryon 
number lies within the interval N m { n < N < N max , where for our model EoS 
the endpoints of iV C rit(^) are N min = 1.49 N Q and N max = 1.78 Nq. As it has 
been shown in jl]], the braking index n{Q) changes its value from n(f2) > 3 
in the region (hi) to n{Q) < 3 in (iv). This is the braking index signal for a 
deconfinement transition introduced by 0]. 




Fig. 4. Braking index due to dipole radiation from fastly rotating isolated pulsars 
as a function of the angular velocity. The minima of n(£2) indicate the appearance/ 
disappearance of quark matter cores. 



In Fig. [| we display the result for the braking index n(Q) for a set of con- 
figurations with fixed total baryon numbers ranging from Nb — 1-55 Nq up to 
Nb = 1.9 Nq, the region where during the spin-down evolution a quark matter 
core could occur for our model EOS. We observe that only for configurations 
within the interval of total baryon numbers 1.4 < Nb/Nq < 1.9 a quark mat- 
ter core occurs during the spin-down as a consequence of the increasing central 
density, see also Fig. 0, and the braking index shows variations. The critical 
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angular velocity ^? C rit(-^s) for the appearance of a quark matter core can be 
found from the minimum of the braking index n(Q). As can be seen from Fig. 
|], all configurations with a quark matter core have braking indices n(Q) < 3 
and braking indices significantly larger than 3 can be considered as precursors 
of the dcconfinement transition. The magnitude of the jump in n(Q) during the 
transition to the quark core regime is a measure for the size of the quark core. 
It would be even sufficient to observe the maximum of the braking index n max 
in order to infer not only the onset of deconfinement (i7 max ) but also the size 
of the quark core to be developed during further spin-down from the maximum 
deviation Sn = n max — 3 of the braking index. For the model EOS we used a 
significant enhancement of the braking index does only occur for pulsars with 
periods P < 1.5 ms (corresponding to fl > 4 kHz) which have not yet been ob- 
served in nature. Thus the signal seems to be a weak one for most of the possible 
candidate pulsars. However, this statement is model dependent since, e.g., for 
the model EOS used by ||, which includes the strangeness degree of freedom, 
a more dramatic signal at lower spin frequencies has been reported. Therefore, 
a more complete investigation of the braking index for a set of realistic EOS 
should be performed. 

3.2 Scenarios with mass accretion 

All other possible trajectories correspond to processes with variable baryon num- 
ber (accretion). In the phase of hadronic stars, fi first decreases as long as the 
moment of inertia monotonously increases with N . When passing the critical 
line N C rit{f2) for the deconfinement transition, the moment of inertia starts de- 
creasing and the internal torque term K mi changes sign. This leads to a narrow 
dip for f2(N) in the vicinity of this line. As a result, the phase diagram gets 
populated for N ~ N clit (Q) and depopulated for N ~ N clit (f2) up to the sec- 
ond maximum of I(N, Q) close to the black-hole line N max (S7). The resulting 
population clustering of compact stars at the deconfinement transition line is 
suggested to emerge as a signal for the occurence of stars with quark matter 
cores. In contrast to this scenario, in the case without a deconfinement transi- 
tion, the moment of inertia could at best saturate before the transition to the 
black hole region and consequently Q would also saturate. This would entail a 
smooth population of the phase diagram without a pronounced structure. 

The clearest scenario could be the evolution along lines of constant fl in the 
phase diagram. These trajectories are associated with processes where the exter- 
nal and internal torques are balanced. A situation like this has been described, 
e.g. by fl30| for accreting binaries emitting gravitational waves. 

In the following we would like to explore which influence the magnitude of 
the external torque K cxt has on the pronouncedness of the quark matter signal. 
In Fig. we show evolutionary tracks (dotted) of configurations in the phase 
diagram of Fig. [I] for different parameter values of the accretion torque and 
different initial values Jo of the angular momentum. 

As we have discussed above, the narrow dip for Q as a quark core signal 
occurs when configurations cross the critical line during a spin-down phase. We 




Fig. 5. Phase diagram for compact stars in the angular velocity - baryon number 
plane with a dividing line for quark core configurations. Trajectories of spin evolution 
are given for different parameter values of the accretion torque dJ/dN in units of [Mq 
km 2 kHz/N©] and different initial values Jo of the angular momentum. The curves in 
the upper left panel correspond to Jo [Mo km 2 kHz] = 300, 400, . . . , 1400 from bottom 
to top. 



can quantify this criterion by introducing a minimal frequency i? m i n above which 
spin-down occurs. It can be found as a solution of the equation for the torque 
balance at the phase border 

dJ/dN = K int (N clit (Q mia ), f2 min )/N . (6) 

The dependence of Sl m i n on dJ/dN shown in Fig. |] can be used to sample 
accreting compact objects from the region in which the suggested quark matter 
signal should be most pronounced before making a population statistics. 

The ideal candidates for such a search program are LMXBs for which the dis- 
covery of strong and remarkably coherent high-frequency QPOs with the Rossi 
X-ray Timing Explorer has provided new information about the masses and ro- 
tation frequencies of the central compact object p0|JTl| |. As a strategy for the 
quark matter search in compact stars one should perform a population statis- 
tics among those LMXBs exhibiting the QPO phenomenon which have a small 
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Fig. 6. Dividing line J? m i n for the deconfinement signal in binary systems with spin-up. 
At a given rate of change of the total angular momentum d J/ dN in units of [Mq km 2 
kHz/No] all configurations with spin frequency £2 m in < O < J? ma x have quark matter 
cores and the population clustering signal is most pronounced. The dashed line is the 
maximum spin frequency J? ma x < Ok{N) for which quark core configurations exist. 

dJ/dN and a sufficiently large angular velocity, see Fig. ||. If, e.g., the recently 
discussed period clustering for Atoll- and Z-sources will correspond to ob- 

jects in a narrow region of masses, this could be interpreted as a signal for the 
deconfinement transition to be associated with a fragment of the critical line in 
the phase diagram for rotating compact stars ||. 

4 Signal for deconfinement in LMXBs 
4.1 Spin- up trajectories for accretors 

We consider the spin evolution of a compact star under mass accretion from 
a low-mass companion star as a sequence of stationary states of configurations 
(points) in the phase diagram spanned by Q and N . The process is governed by 
the change in angular momentum of the star 



-(i(N,n) n) = K eyA , 



(7) 
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where 

K cxt = ^GMM*r - N ont (8) 

is the external torque due to both the specific angular momentum transferee! 
by the accreting plasma and the magnetic plus viscous stress given by -/V ou t = 
K/j 2 r~ 3 , k = 1/3 (3l). For a star with radius R and magnetic field strength 

B, the magnetic moment is given by /i = R? B. The co-rotating radius r c — 

1/3 

(GA//i? 2 ) is very large (r c 3> ro) for slow rotators. The inner radius of the 
accretion disc is 

R , fj,< fi c 
0.52 ta , M > Mc 

where \x c is that value of the magnetic moment of the star for which the disc 
would touch the star surface. The characteristic Alfven radius for spherical ac- 

/ • \ -1 / 7 

cretion with the rate M = raN is ta — I 2/j ~ GMM ) . Since we are in- 
terested in the case of fast rotation for which the spin-up torque due to the 
accreting plasma in Eq. (|§|) is partly compensated by N out , eventually leading 
to a saturation of the spin-up, we neglect the spin-up torque in iV ou t which can 
be important only for slow rotators |]32f , 

From Eqs. (0), (|J) one can obtain the first order differential equation for the 
evolution of angular velocity 

dn_ = K cxt {N, f2)-K iat (N,f2) 

dt I(N,f2) + f2(dI(N,f2)/df2) N ' U 

where 

K int (N, Q) = QN{dI{N, Q)/dN) Q ■ (10) 

Solutions of @ are trajectories in the fl — N plane describing the spin evolu- 
tion of accreting compact stars, see Fig. 0. Since I(N, fl) exhibits characteristic 
functional dependences [^5) at the deconfinement phase transition line A^ cr i t (J7) 
we expect observable consequences in the P — P plane when this line is crossed. 

In our model calculations we assume that both the mass accretion and the 
angular momentum transfer processes are slow enough to justify the assumption 
of quasistationary rigid rotation without convection. The moment of inertia of 
the rotating star can be defined as I(N, Q) — J(N, !?)/!? , where J(N, fi) is the 
angular momentum of the star. For a more detailed description of the method 
and analytic results we refer to Q| and the works of |24|,|25|], as well as p^j27| |. 

The time dependence of the baryon number for the constant accreting rate 
N is given by 

N(t) = iV(to) + {t - t )N . (11) 
For the magnetic field of the accretors we consider the exponential decay |3^] 

B(t) = [B{0) - B^} exp(-t/T B ) + B m . (12) 

We solve the equation for the spin-up evolution @ of the accreting star for 
decay times 10 7 < rs[yr] < 10 9 and initial magnetic fields in the range 0.2 < 
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5(0) [TG] < 4.0. The remnant magnetic field is chosen to be B m = lO^TGQ 
& 

At high rotation frequency, both the angular momentum transfer from ac- 
creting matter and the influence of magnetic fields can be small enough to let the 
evolution of angular velocity be determined by the dependence of the moment 
of inertia on the baryon number, i.e. on the total mass. This case is similar to 
the one with negligible magnetic field considered in ]^. p^.^5| , p6[ where /i < /i c in 
Eq. (H), so that only the so called internal torque term (p_6|) remains. 




N [N J Log(t) [yr] 

Fig. 7. Spin evolution of an accreting compact star for different decay times of the 
magnetic field and different accretion rates. Upper panels: initial configuration with 
7V(0) = 1.4 iV©; Lower panels: JV(0) = 1.55 iV©; J2(0) = 1 Hz in both cases. The 
numbers in the legend box stand for (log(rs[yr]), log(JV[iV0/yr]). For instance (9,-8) 
denotes r B = 10 9 yr and N = W~ a N Q / yr. 

In Fig. Q we show evolutionary tracks of accretors in phase diagrams (left 
panels) and show the corresponding spin evolution Q(t) (right panels). In the 
lower panels, the paths for possible spin-up evolution are shown for accretor 
models initially having a quark matter core (N(0) = 1.55 N®, f2(0) = 1 Hz). 
The upper panels show evolution of a hybrid star without a quark matter core 
in the initial state (N(0) = 1.4 Nq, J?(0) = 1 Hz), containing quarks only in 
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mixed phase. We assume a value of N corresponding to observations made on 
LMXBs, which are divided into Z sources with N ~ 10 -8 -/V©/yr and A(toll) 
sources accreting at rates N ~ 10~ 10 A Q /yr 1 1 3^,|37|. 

For the case of a small magnetic field decay time tb = 10 7 yr (solid and 
dotted lines in Fig.f7|) the spin-up evolution of the star cannot be stopped by the 
magnetic braking term so that the maximal frequency consistent fi max (N) with 
stationary rotation can be reached regardless whether the star did initially have 
a pure quark matter core or not. 

For long lived magnetic fields [tb = 10 9 yr, dashed and dot-dashed lines in 
Fig. 0) the spin-up evolution deviates from the monotonous behaviour of the 
previous case and shows a tendency to saturate. At a high accretion rate (dot- 
dashed lines) the mass load onto the star can be sufficient to transform it to 
a black hole before the maximum frequency could be reached whereas at low 
accretion (dashed lines) the star spins up to the Kepler frequency limit. 



4.2 Waiting time and population clustering 

The question arises whether there is any characteristic feature in the spin evo- 
lution which distinguishes trajectories that traverse the critical phase transition 
line from those remaining within the initial phase. 

For an accretion rate as high as N — 1CP 8 Nq/ yr the evolution of the spin 
frequency in Fig. [?] shows a plateau where the angular velocity remains within 
the narrow region between 2.1 < J?[kHz] < 2.3 for the decay time tb = 10 9 yr 
and between 0.4 < J? [kHz] < 0.5 when tb = 10 7 yr. This plateau occurs for 
stars evolving into the QCS region (upper panels) as well as for stars remaining 
within the QCS region (lower panels). This saturation of spin frequencies is 
mainly related to the compensation of spin-up and spin-down torques at a level 
determined by the strength of the magnetic field. In order to perform a more 
quantitative discussion of possible signals of the deconfinement phase transition 
we present in Fig. ^| trajectories of the spin-up evolution in the P — P plane for 
stars with N(0) = 1.4 Nq and Q(0) — 1 Hz in the initial state; the four sets of 
accretion rates and magnetic field decay times coincide with those in Fig. 0. 

When we compare the results for the above hybrid star model (solid lines) 
with those of a hadronic star model (quark matter part of the hybrid model 
omitted; dotted lines) we observe that only in the case of high accretion rate 
(N = 10~ 8 A Q /yr, e.g. for Z sources) and long-lived magnetic field (tb — 10 9 yr) 
there is a significant difference in the behaviour of the period derivatives. The 
evolution of a star with deconfinement phase transition shows a dip in the period 
derivative in a narrow region of spin periods. This feature corresponds to a 
plateau in the spin evolution which can be quantified by the waiting time r = 

P/P = Q/Q. In Fig. |p (lower and middle panels) we present this waiting time 

in dependence on the rotation frequencies v = Q/(2tt) for the relevant case 
labeled (9,-8) in Figs. @,|. The comparison of the trajectory for a hybrid star 
surviving the phase trasition during the evolution (solid line) with those of a 
star evolving within the hadronic and the QCS domains (dotted line and dashed 
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lines, respectively), demonstrates that an enhancement of the waiting time in 
a narrow region of frequencies is a characteristic indicator for a deconfincmcnt 
transition in the accreting compact star. 



-12 r 




P [ms] 



Fig. 8. Spin-up evolution of accreting compact stars in the P — P diagram. For labels 
and initial values see to Fig. ^. The region of the vertical dashed lines corresponds 
to the clustering of periods observed for LMXBs with QPOs. A dip (waiting point) 
occurs at the deconfinement transition for parameters which correspond to Z sources 
(N = lO _8 A^0/yr) with slow magnetic field decay (tb = 10 9 yr). 

The position of this peak in the waiting time depends on the initial value 
of magnetic field, see Fig. ^[ In the middle and lower panels of that Figure, we 
show the waiting time distribution for B(0) = 0.75 TG and B(0) = 0.82 TG, 
respectively. Maxima of the waiting time in a certain frequency region have the 
consequence that the probability to observe objects there is increased (population 
clustering) . In the upper panel of this Figure the spin frequencies for observed 
Z sources in LMXBs with QPOs JTy] are shown for comparison. In order to 
interprete the clustering of objects in the frequency interval 225 < v\Hz] < 375 
as a phenomenon related to the increase in the waiting time, we have to chose 
initial magnetic field values in the range 1.0 > 5(0) [TG] > 0.6 for the scenario 
labeled (9,-8), see also the dashed lines in Fig|| 

The results of the previous section show that the waiting time for accreting 
stars along their evolution trajectory is larger in a hadronic configuration than in 
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a QCS, after a time scale when the mass load onto the star becomes significant. 
This suggests that if a hadronic star enters the QCS region, its spin evolution 
gets enhanced thus depopulating the higher frequency branch of its trajectory 
in the i? — N plane. 

In Fig. [To] we show contours of waiting time regions in the phase diagram. 
The initial baryon number is iV(0) = 1AN Q and the initial magnetic field is 
taken from the interval 0.2 < B(0)[TG] < 4.0 . 
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Fig. 9. Upper panel: frequency interval for observed Z source LMXBs |_0J; middle 
panel: waiting times r = P/P for scenario (9,-8) and initial magnetic field 5(0) = 
0.75 TG; lower panel: same as middle panel for B(0) = 0.82 TG. Spin evolution of 
a hybrid stars (solid lines) shows a peak in the waiting time characteristic for the 
deconfinement transition. Hadronic stars (dotted lines) and QCSs (dashed lines) have 
no such structure. 



The region of longest waiting times is located in a narrow branch around 
the phase transition border and does not depend on the evolution scenario after 
the passage of the border, when the depopulation occurs and the probability 
to find an accreting compact star is reduced. Another smaller increase of the 
waiting time and thus a population clustering could occur in a region where the 
accretor is already a QCS. For an estimate of a population statistics we show the 
region of evolutionary tracks when the values of initial magnetic field are within 
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0.6 < £?(0)[TG] < 1.0 as suggested by the observation of frequency clustering in 
the narrow interval 375 > i/[Hz] > 225, see Fig. |^. 

As a strategy of search for QCSs we suggest to select from the LMXBs 
exhibiting the QPO phenomenon those accreting close to the Eddington limit 
[ |33"| and to determine simultaneously the spin frequency and the mass p8j for 
sufficiently many of these objects. The emerging statistics of accreting compact 
stars should then exhibit the population clustering shown in Fig. |l0| when a 
deconfinement transition is possible. If a structureless distribution of objects in 
the fi — N plane will be observed, then no firm conclusion about quark core 
formation in compact stars can be made. 
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Fig. 10. Regions of waiting times in the phase diagram for compact hybrid stars for 
the (9,-8) scenario. For an estimate of a population statistics we show the region of 
evolutionary tracks when the interval of initial magnetic field values is restricted to 
0.6 < 5(0) [TG] < 1.0. Note that the probability of finding a compact star in the phase 
diagram is enhanced in the vicinity of the critical line for the deconfinement phase 
transition A r C rit(fi) by at least a factor of two relative to all other regions in the phase 
diagram. 



For the model equation of state on which the results of our present work are 
based, we expect a baryon number clustering rather than a frequency clustering 
to be a signal of the deconfinement transition in the compact stars of LMXBs. 
The model independent result of our study is that a population clustering in 
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the phase diagram for accreting compact stars shall measure the critical line 
N CT it(f2) which separates hadronic stars from QCSs where the shape of this curve 
can discriminate between different models of the nuclear EoS at high densities. 

5 Summary and Outlook 

On the example of the deconfinement transition from hadronic to quark matter 
we have demonstrated that the rotational frequency of accreting neutron stars is 
sensitive to changes of their inner structure. Overclustering of the population of 
Z-sources of LMXBs in the frequency-mass plane {Q — N plane) is suggested as 
a direct measurement of the critical line for the deconfinement phase transition 
since it is correlated to a maximum in the moment of inertia of the compact 
star. 

A generalization of the phase diagram method to other thermodynamical ob- 
servables like thermal and electromagnetic processes is possible and will provide 
a systematical tool for the study of observable consequences of the deconfine- 
ment phase transition. A population statistics in the space of the appropriate 
thermodynamic degrees of freedom allows a direct measurement of the shape 
of hypersurfaces representing the phase border which is characteristic for the 
equation of state of superdense stellar matter. 

In the present contribution we have considered the quasistationary evolu- 
tion of accreting objects where the timescale for accreting about one solar mass 
and thus completing the deconfinement phase transition processes is set by the 
Eddington rate to be about 10 8 yr. Therefore this analysis concerns rather old 
objects. For the cooling process, however, a characteristic timescale is the begin 
of the photon cooling era determined mainly by the relation of neutrino and pho- 
ton emission rates. It occurs typically after 10 2 — 10 4 yr, i.e. for much younger 
objects. Phase transition effects, changing the emission rates, modify cooling 
curves for the neutron star surface temperature typically in this time interval, 
see Voskresensky's contribution to this book j3^| . 

We have focussed here on the deconfinement phase transition only. The phase 
diagram method for developing strategies to investigate changes in the neutron 
star interiors is more general and can be applied to other phase transitions too, 
as e.g. the transition to superfluid nuclear matter f40j , a kaon condensate [|d| 
or color superconducting quark matter p^]43| , f44|| with typical consequences for 
cooling and magnetic field evolution of the neutron star. 

In order to consider the problem of magnetic field evolution in case of color 
superconductivity and diquark condensation one can use the phase diagram 
method to classify the magneto-hydrodynamical effects in the superfluid and 
superconducting phases, e.g. the creation of vortices in the quark core and their 
influence on the postglitch relaxation processes can be discussed Q . 

One can also modify the phase diagram approach to the case of self bound 
rotating quark stars (strange stars) [||). The application of our method could 
visualize the evolution of quark stars from normal to strange stars. 
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The present approach opens new perspectives for the search for an under- 
standing of the occurence and properties of superdense stellar matter in neutron 
star interiors. 
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